Strong suppression of Coulomb corrections to the cross section of e + e pair 
production in ultrarelativistic nuclear collisions 
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The Coulomb corrections to the cross section of e + e~ pair production in ultrarelativistic nuclear 
collisions are calculated in the next-to-leading approximation with respect to the parameter L = 
ln7A7_B (ia,b are the Lorentz factors of colliding nuclei). We found considerable reduction of the 
Coulomb corrections even for large ja^b due to the suppression of the production of e + e~ pair with 
the total energy of the order of a few electron masses in the rest frame of one of the nuclei. Our 
result explains why the deviation from the Born result were not observed in the experiment at SPS 
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Electron-positron pair production in ultrarelativistic 
nuclear collisions is investigated intensively during almost 
two last decades, see recent reviews [3j, This process is 
important in the problem of beam lifetime and luminoc- 
ity of hadron colliders. It is also a serious background for 
many experiments because of its large cross section. For 
heavy nuclei, the effect of higher-order terms (Coulomb 
corrections) of the perturbation theory with respect to 
the parameters Zau and ZbOl can be very important 
(Za and Zb are the charge numbers of the nuclei A and 
B, a w 1/137 is the fine-structure constant). However, 
no evidence of the Coulomb corrections has been found 
in the experiments [l], Q • This circumstance stimulated 
considerable theoretical interest to this process. In the 
set of theoretical works [1, @, 0] it was found that the ex- 
act in Za,bol cross section coincides with that obtained 
in the Born approximation in the ultrarelativistic limit. 
This statement was considered as an explanation of the 
experimental results [HQ- However, this conclusion con- 
tradicted to the result obtained in Ref. [§[ with the help 
of the Weizsacker- Williams approximation in the lead- 
ing logarithmic approximation. This contradiction has 
been resolved in Ref. [9(. It was shown that the wrong 
conclusion of Refs. [lHI 0] on the absence of Coulomb 
corrections was due to the bad treatment of condition- 
ally convergent integrals. Consistent approach of Ref. [9| 
results in the Coulomb corrections which coincide with 
those from Ref. @ . Thus absence of the Coulomb correc- 
tions in the experiments [J 0] has remained unexplained. 
In the present paper, we explain this experimental result 
by the suppression of the Coulomb corrections due to the 
account of the next-to-leading term. 



W{b) = np n(b) , °t = / d 2 bW(b) = na n , 

1 ** 1 



where a n = j d 2 bP n (b) is the cross section of n-pair pro- 
duction. The cross section ot can be presented in the 
form: 



ab 



(2) 



where er° oc (Z aol) 2 (Z bol) 2 is the Born cross section, 
a A and <r B are the Coulomb corrections with respect 
to nucleus A and B, respectively (containing the terms 
proportional to (ZbOl) 2 {Zao) 2ti and {Z B Q) 2n {Z aol) 2 , 
n 5* 2), and a AB is the Coulomb corrections with re- 
spect to both nuclei (containing the terms proportional 
to {Z Ba) n (Z acx) 1 with n, I > 2). The cross section a 
coincides with the Born cross section of one pair pro- 
duction, which was calculated many years ago in Refs. 
[Tol . [TlT | . In the leading logarithmic approximation, the 
quantities a A ' B oc L 2 and a AB oc L were obtained in 
Refs. and Ref. [HI], respectively. 

The leading logarithmic approximation for W(b) pro- 
vides the factorization of P n (b) 0QHE1, 

so that 



Pn(b) = 



W n (b) 



-W(b) 



(3) 



The function W(b) was calculated in Refs. [17], Gi 



2C . 12111 in the Born approximation and in Refs. 
23,r24L 251 with the Coulomb corrections taken into 



account. Using Eq. ([3]), the cross section u\ of one pair 
production can be represented as a sum of <jt and the 
unitarity correction <r un i t 



Since the nuclear mass is large compared to the elec- 
tron mass, it is possible to treat the nuclei as sources 
of the external field and calculate the probability P n (b) 
of n-pair production at a fixed impact parameter b. It 
is convenient to introduce the average number W(b) of 
produced pairs and the number-weighted cross section 
or as 



Cunit 



(Tl = a T + Cunit , 
d 2 



bW(b) ( 



1-e 



-W(b) 



(4) 



The existence of the unitarity correction was first recog- 
nized in Ref. [26[ (see also review Q). Numerical evalu- 
ation of this correction was performed in Refs. [20L [27l|. 
The main contribution to a\ is given by the term a in 



2 



ctt, Eq. ([2]), and is known with high accuracy [ToL [Til] . 
The terms a A and a B in also give important contri- 
butions to a± . In the leading logarithmic approximation, 
these terms have been derived in Refs. [8, 9]. The last 
two contributions, a AB and a un a, to a\ are rather small, 
see Refs. [Hill. 

In the present paper, we calculate the leading cor- 
rections to a A ' B (which are also the corrections to o\). 
We show that these corrections essentially diminish the 
magnitude of a A ' B even for the parameters of LHC 
(lA — 1b ~ 3000, Za = Zb= 82). It is convenient to 
calculate a A in the rest frame of the nucleus A, where the 
nucleus B has the Lorenz factor 7 = 27,473 at ^a.b 3> 1- 
Note that a A , being proportional to (Zb&) 2 , can be di- 
rectly calculated as the Coulomb corrections to <Ti with 
respect to the parameter Z A a, so that it can be repre- 
sented as 



duo 



dQ 2 



dujdQ 2 



dnj_(uo, Q 2 ) 
diodQ 2 



a ± (oo,Q 2 ) 



-a {l (uo,Q 2 ) 



(5) 



Here 



dn±{uo, Q 2 ) 



dnn(uj,Q 2 ) 



1 



(oo/j) 2 \ diodQ 2 



7T V Q 2 

Z%a du dQ 2 



00 Q 2 



7T 



W Q 2 



(6) 



are the numbers of virtual photons 7^ 1, with the en- 
ergy uo, the virtuality — Q 2 < 0, and the transverse and 
longitudinal polarizations, respectively. The quantities 
a±(uo, Q 2 ) and a\\(u>, Q 2 ) are the Coulomb corrections to 
the cross sections of the processes 7^ «A — > e + e~ A. 

Let us discuss the contributions to a A of different re- 
gions of integration with respect to uo and Q 2 . 

The leading logarithmic contribution oc L 2 comes from 
the integration of aj_ over the region 



I. m<^uo<^mj, (w/7) 2 < 



f < m 2 . 



(7) 



The leading correction oc L comes from the following 
regions: 



II. 
III. 
IV. 



Q 2 
Q 2 



m , ni uo 7m 

\2 



m<u< 7m 



m, (m/7) 2 <Q 2 < to 2 



(8) 
(9) 
(10) 



Note that the cross section mi gives logarithmically en- 
hanced contribution only in region II. Therefore, since 
we are going to keep the terms proportional to L 2 or L, 



we can write a A as 
a A = at + 5a A 



dQ 2 



dn±(uo, Q 2 ) 



diodQ 2 



o-j_(oo, Q 2 



dn\\ (uj, 



diodQ 



5a 1 



, a||(oo,Q 2 ) 

OO 

2 dn±(cv, Q 2 ) 



(ii) 



dQ 2 



2m 



diodQ 2 



S<j ± (lo,Q 2 ) 



5a ± (uo, Q 2 ) = a±(uo, Q 2 ) - er±(oo, Q 2 ) . 



(12) 



The quantities erjii (00, Q 2 ) can be calculated within 
the quasiclassical approximation. Following the method 
described in detail in Ref. 1231, we obtain 



(00, Q 2 



a 



Re 



2e- 



P2 



ke) D_ 



de Sp 
2e* • 



Pi 



dr l( 2r 2 e- lkr 
- ke*) D 



D_ = D(r 2 ,ri|e), D+ = D(r x , r 2 | £ - w) , (13) 

where _D(r2,ri|e) is the quasiclassical Green function of 
the squared Dirac equation, e = (0,1,0,0) for a± and 
e = (0, 0, 0,Q/uj) for an in the frame where k is directed 
along z axis. Using the explicit expressions for the Green 
functions from Ref. (28[, we obtain the results for these 
cross sections: 



ct±(oo,Q 2 ) 



n.V / dy 





1 + 2(1- 2yy) (1 + yyQ 2 /m 2 



a«(oo,Q 2 ) = AaN / dy 



(1 + yyQ 2 /m 2 ) 2 
y 2 y 2 Q 2 /m 2 



(1 + VvQ 2 /m 2 ) 



2\2 



N 



y = 



MZ A af 
3m 2 

i-y, 



(14) 



where tp(x) = d\nT(x)/dx. These formulas agree with 
the result of Ref. [2i| if one takes into account the missing 
factor yy in an pointed out in Ref. [30j . Substituting Eq. 
(TTJJ in Eq. (Ti"Tj) and taking the integrals over w and Q 2 , 
we obtain within the logarithmic accuracy 



7(Z B a) 2 N 
3tt 



20 
— j 
21 



(15) 



We remind that L = ln(j A j B ) = 111(7/2). The result (JUJ 
is in agreement with those obtained in Refs. [30l [3l| . 

Let us pass to the contribution Sa A , Eq. (|12p. which 
was not considered so far. In Ref. '31] it was conjec- 
tured that the term Sa A can be safely omitted. We show 
below that this guess is completely wrong. The func- 
tion 5a ± (00, Q 2 ) in the integrand provides convergence 
of the integral over to in the region to ~ m. The log- 
arithmically enhanced contribution is given by the re- 
gion (m/7) 2 <C Q 2 -C m 2 of integration over Q 2 . Since 
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FIG. 1: The dependence of G(Za) on Z. 



The function G(Za) is shown in Fig. [TJ It is seen 
that G(Za) varies slowly from —6.6 for Z = 1 to —6.14 
for Z = 100 being large for all interesting values of Z. 
The large value of G leads to a big difference between 
a A from Eq. (|17[) and its leading logarithmic approxi- 
mation a A A = 7 {Z sot) 2 N L 2 /(3ir) even for very large 7. 
This statement is illustrated in Fig. [2] where the ratio 
<tj4 / <t la i s snown as a function of 7 (solid curve). If one 
omits the contribution 8a A and use a A s , Eq. ([15]) . as an 
approximation to <j A , then the contribution of linear in L 
term becomes much less important, see the dashed curve 
in Fig. O Note that for Pb-Pb collisions at LHC one has 



1.2 — - 



Q 2 <C rn 2 in this region, we can make the substitution 
<5ct ± (w,Q 2 ) -> <5ct ± (w,0) in Eq. Then we take the 

integral over Q 2 and obtain 



x A 
G(Z A a) = 2 



7(Z B a) 2 NG{Z A a) 



3tt 



dui 

2m U 



<7j_ (OJ, 0) 
<7± (00,0) 



(16) 



The quantity a±(u>,0) = ct^a(u)) is the Coulomb correc- 
tions to the cross section of e + e~ pair production by real 
photon in the Coulomb field, and a± (00, 0) = 7aN/3. 
Taking the sum of Eqs. (Tl"5|) and (TIT)]) , we finally obtain 



a in the next-to-leading approximation 



7(Z B a) 2 N 
3^ 



L + G(Z A a) 



20 
21 



L 



(17) 



In order to calculate the function G(Z A a) it is neces- 
sary to know the magnitude of the Coulomb corrections 
g iA (uj) in the energy region where the produced e + e~ 
pair is not ultrarelativistic. The formal expression for 
it, exact in Z A a and ui, was derived in Ref. 32]. This 
expression has a very complicated form causing severe 
difficulties in computations. The difficulties grow as u> 
increases, so that numerical results in Refs. [32l. l33j were 
obtained only for u> < 5MeV. In a set of later publica- 
tions [H [H, [H, H3 (see also reviews [MUl) the mag- 
nitude of a lA (uj) has been obtained for higher values of 
lo and several Z A . In the high-energy region u> 3> m, the 
consideration is greatly simplified. As a result, a rather 
sim ple form of the Coulomb corrections was obtained in 
[40L l4l[| in the leading approximation with respect to m/oj 
and in Ref. [2^] in the next- to- leading approximation. In 
Ref. [42| . a simple formula, which correctly reproduces 
the low-energy results and the high-energy limit, was sug- 
gested. This "bridging" expression has high accuracy at 
intermediate energies and differs from the exact result for 
a lA (uj) only in the region close to the threshold lu = 2m. 
For our purpose, this difference is not important because 
in this region the ratio a lA {uj) / a lA (po) can be neglected 
in comparison with unity. 
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FIG. 2: The ratio a A /a A A (solid curve) as a function of 7 for 
Z A = 82. Here a£ A = 7(Z B a) 2 NL 2 /(3tt) is the Coulomb cor- 
rections calculated in the leading logarithmic approximation. 
Dashed curve shows the ratio & A /a A A . 



7 w 1.8 X 10 7 and o A ja A K w 0.66. ForAu-Au collisions 
at RHIC one has 7 w 2.3 x 10 4 and cr A /a A A w 0.42. 
For the experiments at SPS [l], the Lorentz factor 
was 7 w 200. Naturally, we can not use the result (fl7|) 
obtained in the logarithmic approximation in the region 
7 < 500 where the logarithmic correction to a A becomes 
larger than the leading term <J A A . However, we can claim 
that, due to the strong compensation between the leading 
term and the correction, the Coulomb corrections a A are 
much smaller than a A A at 7 < 500. Therefore, this natu- 
rally explains why there was no evidence of the Coulomb 
corrections in the experiments P, Q • 

Let us discuss now the importance of the Coulomb cor- 
rections (7 in comparison with the Born cross section 
cr°. The ratio a A /(T° is shown in Fig. [31 In the next- 
to-leading approximation for a A this ratio (solid curve) 
is small (< 5%), while the same ratio obtained with a A 
approximated by a A A reaches 20% at 7 ~ 1000. 

To summarize, we have calculated the Coulomb cor- 
rections a A to e + e~ pair production in the next-to- 
lcading logarithmic approximation. After the account 
of the next-to-leading term, the magnitude of a A be- 
comes small in comparison with the Born cross section, 
in contrast to the leading term cr A A . The big difference 
between our result and previously suggested one has a 
simple explanation. The latter was based on the use of 
the high-energy asymptotics for the Coulomb corrections 
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FIG. 3: The ratio a A /a° (solid curve) as a function of 7 for 
Z a = 82. Dashed curve shows the ratio cr^/cr . 



to the photoproduction cross section instead of the exact 
Coulomb corrections. However, the exact Coulomb cor- 
rections are strongly suppressed in a rather wide region 
2m < u> < 20m. Note that our results, combined with 
a AB , Ref. [12], complete the calculation of linear in L 
terms in the number- weighted cross section <jt- 
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